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Abstract 

Within the framework of Gaussian equivalent representation method a new procedure of ob- 
taining equations of state for simple liquids is discussed in some technical details. The developed 
approach permits one to compute partition and distribution functions for simple liquids with arbi- 
trary form of the central two-body potential of inter-molecular interaction. The proposed approach 
might become of great use for computing thermodynamic and structural quantities of simple parti- 
cle and polymer systems. We believe that this technique can also provide an interesting possibility 
to reduce the sign problem of other methods of computer simulation based on a functional integral 
approach. 
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I. INTRODUCTION 



The calculation of partition and distribution functions is a basic problem of statistical 
pliysics[l|. All thermodynamics characteristics of statistical systems are determined by these 
functions. As is well known, the calculation of those quantities is a formidable problem^, 

aa. 

In this paper we develop a functional integration method for systematic approximate 
calculations of classical partition functions of two-body potentials with positive and negative 
Fourier transforms over the entire density and temperature range. 

The Gaussian equivalent representation method has been recently introduced by Efimov 
and Ganbold in the context of quantum-field theory and statistical physics to compute 
integrals over Gaussian measure y, 0]. The GER approach has already been proven to be 
very effective for computing thermodynamic properties and structural quantities of simple 
classical many-particle systems interacting with purely repulsive potentials like the Gauss- 
core or the Yukawa potential, possessing positive Fourier coefficients jsl, [oL 10|as well as for 

12| . Moreover, 



calculating the thermodynamics properties of flexible polymer systems [jjj 



it has successfully been ernp 



Monte Carlo simulation 
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oyed to reduce the numerical sign problem in conjunction with 



14| . In the meantime, the real liquid description demands 



3,Q- 



of considering the potentials having both attraction and repulsion parts IJJ In the 

present work the author extends the approach for systems, where the particles interact 
through potentials with positive and negative Fourier coefficients. This increases the range 



of applicability of this method for computational simulations 



19 
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221. 



II. THE EQUATIONS OF STATE IN THE THEORY OF SIMPLE LIQUIDS 

The simple fluid of the particle density n = N/V is thought of as a dense cloud of 
particles occupying the volume V and interacting via two-body potential of the form 
Thermodynamics of such a system is described by the partition function 

For given two-body potential V{'x — x') the free energy of the system can be computed and 
can be written in the form 



2 



E(n,P) = -4 lim InZy. 

P \V\^oo 



(2) 



For systems of particles interacting via the two-body potentials having attraction Vi and 
repulsion V2 parts the total potential can conveniently be represented in the matrix form 



Vf^ 



V, 



-1 



The differential operator V ^ (x — y ) satisfies the equation 



J dyV-\^ - y)V{y - x') = <5(x - x'), 



and has a form 



V-^(x-y) 



V{-A) 



^(x-y) 



or V = I 



(27r)3 v{k) 



where V(k) — J dyiV{x)e^^^ is the Fourier-image of the potential. The identity 



£^^-^i<pv-^<p)+i{b4>) ^ g-i(bVb) 



with 



TV 



6(x) 



/^^<5(x-Xj) 
permits one to represent Boltzmann factor as 

Cv 



e-/?Er<,-^(xi-xj) ^ givfy(o)-|E5i^(xi-xj 

N 



,iV^</.(xj)+|y(o) 



^ .7=1 



iV;g</>(xj): 



where the following notation is utilized 



(ix6(x)0(x), 
(bVb) =11 rfxrfy6(x)y(x-y)6(y), 

= / / dxdy0(x)y-^(x-y)0(y). 



Let us consider the partition function in the form of functional integral 



Zv = 



V 



dxi 

W\ 



(ix 



(3) 

(4) 
(5) 

(6) 



(7) 

(8) 
(9) 
(10) 



\V 
1 

W 



N 



VdetVi J y/det\^ 



g-i(<AyrV)g-i(V'V2-V) 



1 N 



V 



V 



VdetV 



^_i(<I.y-i$) 



/ dx: 

Jv 



N 



N\ 1 



2m\V\^ T z^+^ 



dz 



Mz) (11) 
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where 



$(a:) = (0(a:),V^(x)), J=( ' I, /y(^) = / ^^e-i(*^-^*)+^/v'^'^^^'^*'-. (12) 



—I 



VdetV 

Following the line of argument of the GER, the integral Iv{z) is written as follows 



'^^^^ y detvj v/d^ ^ ^ 



VdetD 



where 



(JVJ) = J dxj rfyJ5(x - x')l^(x - y) J5(x' - y) = V\(0) - ^2(0). (15) 
As a result, we obtain two equations 

Equation I : i j j dy($(x) [^-^(x - y) - ^-^(x - y)]|<f (y)) 

+ i y t/xy rfy[-/?($J)2]e*^(*0'^)+^['^(^-^)^l = (16) 

or, more specifically this equation reads 

JK. ik(x-y) I- JK. ik(x-y) . A, . 

(27r)3 v^(k) y (27r)3 /^(k) '"^ ^ 

J— = n2e^v/3(1'oJ)+§[J(V-D)J]^ (17) 

V{k) D{k) 



Equation II: -($^-^$0) + i^v^^ rfx($(x) J)e'^(*"^)+2 [^(^-^)-^] (18) 



or 



*„.^.V^f/.x(*.)e.v^<-Hl.--«..^1.2. (19, 

Assume $0i = iCi/P, where 

a = zpVijiO)Jje~''-'^^^^-'^''-''^'\ (20) 



Making use of (!20|) in (fT7|) . we obtain 

1 1 \ Ci / 1 1 \ Ci 



Vik) D{k)J,, Vm Wik) Dik)J,, V\(0) 

1 1 \ / 1 1 \ C2 



V{k) D{k)J,, \V{k) D{k)J,, ^(0) 



(21) 
(22) 
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This can be combined in the matrix 



1 



1 



+ 



C2 



D{k) 

Thus, we obtain for D{k) 



Vi(k) Vi{0) V2(0) 

C2 1 Cl 



D~\k). 



(23) 



V2(0) 



V2(k) Vi(0) 



1 



Cl 



C2 



V2(fc) yi(o) 



^2(0) 



C2 



1 



+ 



Cl 



(24) 



where 



A = det 



D{k) 



+ - 



V2(0) Vi{k) Vi(0) 
Cl \ / 1 Cl 



Vlik) Vi{k)J \V2{k) \/i(0) 



C2 



^2(0) 



At large we can use Stirhng's formula 
N\ 



\V\^ 



g|y|{ninn-n) _ ^\v\h ^ ^^^^^ = nlnn - n. 



It is convenient to introduce function 



\V\R{c) = \V\fn + Wo-N\nz 



(25) 



(26) 



where 
Wo 



1, detD 1 
- m 

2 detV 



\V\ 



<J)2 

^01 



dk 



2 2 Jy 

A-\k) 



In-^ 



<J)2 

^02 



Viik)V2ik) 



+ Tr{l- D{k)V~\k)) 



2Vi{0) 2^2(0) f3V,{0) 



and 



In^ = In + cJ- -(J[V - D]J). 

/3Vi(0) 2^ ^ ^ ' 

After some algebra we obtain the partition function in the form: 

1 



2m J z 



<^^\V\R{c(z)) 



VdeW 



where 



t> Jv 



D- 



(27) 



(28) 



(29) 



(30) 



The partition function integral has been derived in a similar from in section III of reference 



231] for potential models with positive Fourier coefficients (see Eqs. (29)-(31)). As a point 



of interest, we consider the lowest approximation partition function having the form 

1 f dz 



7° 



2m J z 



(31) 



In the case V — > oo integral (13T]) could be done using the steepest descent method 



70 



27[i I c 



The point of global maximum is defined by the equation 

^R{c) = h{c) + hie) + hie) + N, 
ac 

The explicit form of functions /i(c), /2(c), /3(c), reads 



(32) 



(33) 



h{c) 

h{c) 
h{c) 
Nic) 



°° dk 
(2^ 

°° dk 



k' 



1 



1 



Vi{k)V2{k) V Vi{o) 



{2nyVi{0)A^{k)\V2{k) V^ik) 

dk Pk"^ n /I 1 y 



2 + ^ 



^i(A:)-\/2(A;) 



n i - 



PViiO) c ' y Vi{0) J V^\0)(3 \ Vi{0) J 
The final expression for R{c) is given by 



R{c) = Mi{c) + M^ic) + P{c) 



(34) 



where 
Mi(c) 

M2(c) 



(27r)2 
dk 



(Vi{k) - V2{k) 



1 



(27r)^ 



In- 



A{k)V,{k)V2{k) 
1 



A(A;)Vi(A;)\/2(A;) A(A;)Vi(A;)V2(A;) 



2 + ^ 



P(c) = n ( Inn — 1 — In 



+ 



— nc ] + 



c 



All other thermodynamic functions may be found from E{n,j3) by the Maxwell relations in 
thermodynamics 



p__(dE\^ _l^\nZ{n,(3) 



\dV ^ T 

In particular, from above it follows 



dV 



R{n,P) = P{n,P). 



S 



dE 
df 



(35) 



(36) 



where P is pressure of system. 
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III. SUMMARY 



The developed procedure of computing the equation of state can be summarized as fol- 
lows. First, we solve equation ([33D, and, second, the obtained roots are inserted in the 
equation (IMIl which is the equation of state. Thus, the developed procedure permits one to 
get the equation of state for simple liquids, composed of particles interacting via two-body 
potential with attractive and repulsive counterparts and having bound states. The applica- 
tion of this procedure to the simple liquid models with specific potentials is the subject of 
forthcoming article. 

We hope that the techniques presented in this Letter can also be useful for other fields 
of computer simulation, where the sign problem does occur. We can contribute in this way 
to establish the auxiliary field methodology as a standard tool for computation and this 
technique can also provide an interesting possibility to reduce the sign problem of other 
methods of computer simulation based on a functional integral approach. 

This work has been initiated by discussions with Prof. Garry V. Efimov to whom author 
is very indebted. 
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V. APPENDIX 

A. Normal form of functional 

Let us introduce normal form of functional with respect to Gauss measure. Consider the 
equality: 




(37) 



and identity: 




(38) 
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The normal form of functional with respect to Gauss measure with Green function D will 
understand multiplication 

This definition is identity and valid with any b, therefore expanding both hands of Eq. [39] 
with respect to b we obtain 

:0(x):z) = 

:0(xi)0(a;2):D = </)(xi)0(x2) - D{xi,X2), 
:0(a;i)0(x2)0(x3):i) = 0(a;i)0(x2)0(x3) - D{xi, X2)(p{xs) - D{x2, X3)(f){xi) - D{x3, Xi)0(a;2), 



We can use the functional in normal form as: 

da^^D-e'^'^^-D^l (40) 

da^^D'-4>{xi <j){xn)y-D = 



In particular: 



where 



/(^) = J = J rfa^^^e^^l^l (41) 

with normalization: 

m = I §^e-i(^^-V) = I da,,^ = 1. 

Any path integral over Gaussian measure for analytical can be written in normal form 
with some Green function D: 

W[(t)] = J rf/i^,e'(''^) = J rf/i^,e-5(^^'');e*('"^):o (42) 

This path integral convenient to represent in the form: 

W[<j)] = Wo + t{Wi<f)) - ^\(j)W2(l)-.D + '■.Wi[<j)]-.D 



Where 
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el = e''-\-z . 

2 2 



B. Equations 

Let us consider path integral (HTl) and make use of the following equivalent transforma- 
tions. First, we shift the variable of integration </)(x) — > Second, we write the 
functional of interactions in the normal form with respect to Gauss measure for new kernel 

B^^ixx.X'i)^ we obtain 

g ^{^^B-^4>)^-\(^4,{A-^^B-^\4>)-{4,A-^li,)-\{iA-^i)+gW{4,+i\ 



Ca J Cb 



C 



e 



A 



gWo+ig{Wi(P)-^:<t>W2<t>:B+9-WM:B 



(43) 



The major contribution to the functional integral gives Gauss measure da^^B, therefore the 
linear and quadratic terms over the integration variable 0(a;) should be absent. Thus we 
obtain two equations 

Equation I : -{(j)A-^^) + ig{Wi^) = 0, (44) 
Equation II : -^MA'' - 5-^):^ " |K0W^20);b = 0. (45) 



For more details see 



l6|. 
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